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Abstract 


Theoretical study was carried out to investigate the possible improvement of SOFC performance by using a non-uniform potential operation 
(SOFC-NUP) in which the operating voltage was allowed to vary along the cell length. Preliminary results of a simple SOFC-NUP with a cell 
divided into two sections of equal size in term of range of fuel utilization (Uf) indicated that the SOFC-NUP can offer higher power density than 
an SOFC with uniform potential operation (SOFC-UP) without a reduction of the electrical efficiency. In this work, voltages and section splits 
were optimized to obtain the maximum power density of the SOFC-NUP. At the optimum splits (S,; =0.55 and S,=0.45), the power density 
improvement as high as 9.2% could be achieved depending on the level of electrical efficiency. It was further demonstrated that the increase in 
the number of separated section (n) of the cell could increase the achieved maximum power density but the improvement became less pronounced 


after n>3. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Fuel cells are promising electrochemical devices that con- 
vert the chemical energy of a fuel directly into electrical energy. 
Compared to conventional power generation processes, fuel 
cells are attractive due to their better environmental friendli- 
ness, practical noise-free operation, and higher efficiency. A 
number of researches have been carried out in many directions 
with a major aim to improve performances of the fuel cells and 
their systems. For instance, several gas turbine cycles such as 
steam injected gas turbine cycle, gas turbine/steam turbine com- 
bined cycle, and humid air turbine cycle have been combined 
with solid oxide fuel cell (SOFC) for efficiency improvement 
[1]. An operation with anodic offgas recirculation was proposed 
for a polymer electrolyte membrane fuel cell (PEMFC) system 
equipped with a fuel processor. Under this operation, a signif- 
icant efficiency increase for the fuel processor and the gross 
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efficiency of the combined system of 30% were reported [2]. 
A novel gas distributor with current-collecting elements dis- 
tributed in gas-delivery fields for effective current collection and 
heat/mass transfer enhancement was designed to improve power 
density [3,4]. In addition, many researches have been devoted 
on development of cell components with superior characteristics 
[5-7]. The performance improvement of PEMFC from operation 
at a higher pressure was reported although the power loss due to 
air compression was taken into account in the comparisons [8]. 

Among several procedures for improving fuel cell perfor- 
mance, the use of non-uniform potential operation concept for 
fuel cells in which the cell voltage is allowed to vary along the 
cell length is another interesting approach. However, until now 
only some works have focused on the potential benefits of this 
approach. It was demonstrated that an improvement in electri- 
cal efficiency of about 1% could be achieved by splitting the 
cell of a molten carbonate fuel cell (MCFC) into two sections 
[9]. Selimovic and Plasson [10] examined performances of net- 
worked solid oxide fuel cell (SOFC) stacks combined with a gas 
turbine cycle. Two multistage configurations, i.e. (i) both anode 
and cathode flows were serially connected and (11) only the anode 
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Nomenclature 

a constant in Eq. (11) (Q m) 

A area (m2) 

b constant in Eq. (11) (K) 

E activation polarization energy in Eqs. (12) and 
(13) (kJ mol!) 

Eo open circuit voltage (V) 

F Faraday constant (96485.34) (C mol`!) 

i current density (A m~?) 

LHV _ lower heating value of methane feed (W) 

m constant polarization parameters in Eqs. (12) and 
(13) 

n number of separated section 

p partial pressure (atm) 

P power density (W cm~?) 

r activation polarization parameters in Eqs. (12) 
and (13) (A m7?) 

R universal gas constant (8.31447 x 107°) 
(kJ mol~! K7!) 

Sp section split 

T absolute temperature (K) 

Uf fuel utilization (%) 

V operating voltage (V) 

W electrical work (W) 


Greeks letters 


ô thickness (m) 

E electrical effiency (%) 
n overpotential (Q m?) 
Q potential (V) 

p specific ohmic resistance (Q m) 
Subscript 

A anode 

Act activation 

C cathode 

Conc concentration 

k section number 

Ohm ohmic 


flow was serially connected while the cathode flow was parallel 
connected, were considered. An increase of system efficiency 
of about 5% was reported for the former configuration mainly 
due to an improved thermal management. The similar multi- 
stage configurations were also considered for a combined heat 
and power MCFC plant [11]. Detailed flowsheet calculations 
showed that the improvement in efficiency was about 0.6% for 
the former configuration, and 0.8% for the latter configuration. 
The concept was also extended to PEMFCs divided into many 
stages (or stacks) of equal size [12]. It was demonstrated that 
the non-uniform cell potential operation allowed for enhanced 
maximum power densities compared to the traditional concept 
involving a uniform cell potential distribution. The improvement 
within 6.5% was reported. 


In this study, the concept of non-uniform potential operation 
implemented to SOFCs fed by methane was investigated in an 
effort to optimize cell operating voltages and sizes so that a max- 
imum power density could be achieved without a reduction of 
electrical efficiency. Values of power density of a simple SOFC 
with a cell divided into two sections whose operating voltages 
and sizes were allowed to vary were compared to those of a 
typical SOFC with uniform cell potential at various electrical 
efficiencies. Additionally, the effect of number of cell section 
on the obtained performance was determined. 


2. Theory 


The schematic diagram of an SOFC with non-uniform poten- 
tial operation (SOFC-NUP) is illustrated in Fig. 1. Compositions 
of fuel and air streams change along the cell channel according 
to changes in value of fuel utilization defined as the mole of 
hydrogen electrochemically consumed divided by the theoret- 
ical mole of hydrogen generated from complete reforming of 
the methane feed. In section k, the fuel utilization changes from 
Ugy—1 to Ugy, and the cell is operated at a constant potential of 
V within the section. In practice, the non-uniform cell poten- 
tial can be realized by using segmented current collectors for a 
single-cell SOFC or it can be carried out in a series of SOFC 
stacks operated at different stack potentials. The section split of 
section k (Spx) is defined by 


~ Utk — ei 


S k 
P: U 
f, final 


(1) 

When the SOFC is fed by a non-hydrogen fuel (e.g. methane), 
a reformer is generally required to reform the fuel with an oxi- 
dant (e.g. steam) to a hydrogen-rich stream before feeding to 
the SOFC stack. The main reactions involved in the production 
of hydrogen from methane and steam are the methane steam 
reforming and water gas shift reactions as shown in Eqs. (2) and 
(3), respectively. 


CH4 +H20 = 3H2+CO (2) 
CO + H20 = H2 + CO2 (3) 


In the SOFC stack, both hydrogen and CO can react electro- 
chemically with oxygen. However, it is assumed that the CO 
electro-oxidation is neglected. It was estimated earlier that about 
98% of current is produced by H2 oxidation in common situa- 
tions [13]. This is due to the fast rate of water gas shift reaction at 
an SOFC operating temperature. The electrochemical reactions 
of hydrogen and oxygen take place according to Eqs. (4) and 


(5). 
Anode: H +0% = H20 + 2e7 (4) 
Cathode: O2 +4e7 = 20% (5) 


Electromotive force (Eo) of the cell is a difference of potentials 
between both electrodes of the cell. It can be represented as 
follows: 


Eo = |c — gal (6) 
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Fig. 1. Schematic diagrams of SOFC-NUP. 


where gc and ga are the potentials of the cathode and the 
anode, respectively. The electrode potential can be calculated 
from Nernst equation which can be expressed as follows: 


RT l 7 
g= 4F n po, (7) 


where R is the universal gas constant, T the absolute temperature 
and F is the Faraday’s constant. The partial pressure of oxygen 
in the cathode chamber is calculated from its mole fraction while 
the partial pressure of oxygen in the anode chamber is given by: 


2 
PH20 
= — 8 
Po, ( zao) (8) 


where p; is the partial pressure, and K is the equilibrium constant 
of the hydrogen oxidation reaction. 

When the SOFC is operated at a potential of Vz, the local 
current density can be determined according to Eq. (9). It should 
be noted that as the value of Eg changes along the cell length 
within the section due to the change in gas compositions, the 
current density varies within the cell section. 


: Eo — Vx (9) 
— 
NOhm + NAct,A + NAct,C + NConc,A + NConc,C 
where 
Nonm = > pj8j (10) 
pj = aj exp(bjT) (11) 
=j 
4F Po m — Eç 
TE faeces B02 eis 12 
NAct,C ($r p exP | -RF (12) 
2F PH i —EA 7 
EEN (acai peasy eens: 13 
NAct,A (iral z exp | T (13) 


The parameters used in the calculations are provided in Table 1 
[14] and Table 2 [15]. It was reported that the correlations of the 
activation polarization can predict the cell performance close 
to the Butler—Volmer equation within the temperature range of 
1173-1273 K [16]. In this study, the operating temperature is 
kept at 1173 K which is in the reliable limit. 

To simplify the calculations of the SOFC performance, it is 
assumed that both fuel and oxidant are well-diffused through the 
surface of the electrodes. Therefore, concentration polarization 
losses (HConc,A and NConc,c) can be omitted. This assumption 


Table 1 
Resistivity and thickness of cell component [12] 


Material used Ni-YSZ/YSZ/LSM-YSZ 


Anode thickness (wm) 150 

Anode ohmic resistance constant a=0.0000298, b = — 1392 
Cathode thickness (um) 2000 

Cathode ohmic resistance constant a=0.0000811, b=600 
Electrolyte thickness (wm) 40 

Electrolyte ohmic resistance constant a=0.0000294, b = 10350 


Interconnect thickness (wm) 10 
Interconnect ohmic resistance constant a=0.001256, b = 4690 


is valid when the cell is not operated at too high current den- 
sity or too low concentration. In addition, it is further assumed 
that the cell is operated at isothermal condition and the fuel 
stream is always at its equilibrium composition along the length 
of the SOFC cell. It should be noted that an external reformer 
is usually connected to the SOFC system in order to generate 
hydrogen-rich feed for the stack and to suppress the cell deacti- 
vation due to the carbon formation. In addition, state-of-the-art 
SOFC nickel cermet anodes are usually active for the steam 
reforming and shift reactions particularly at high operating tem- 
perature of SOFC [16,17]. Calculations of the thermodynamic 
equilibrium composition are accomplished by following details 
given in our previous work [18]. 

When the current density is known, the cell area (Ax) and 
electrical power (W,) involved in each cell section can be deter- 
mined. The values of overall electrical efficiency (£) and average 
power density (P) can be determined according to Eqs. (14) and 
(15), respectively. 


n W 
k=1 
n 
W 
p = k= Wi (15) 


where LHV is the lower heating value of methane feed. 


Table 2 
Summary of activation polarization parameters [13] 

r (Am?) Ex pol (kJ mol™!) m 
Cathode 1.489 x 1010 160 0.25 
Anode 2.128 x 108 110 0.25 
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3. Results and discussion 


Fig. 2 shows the performances of a typical SOFC with uni- 
form potential operation (SOFC-UP). Each line represents the 
results at a constant value of fuel utilization (Uf). It should 
be noted that because concentration polarization becomes an 
importance loss at high current densities and small fuel concen- 
trations (Uf > 80%) [15], the simulations were performed only 
in the ranges of high voltages (low current densities) and fuel 
utilizations (Uf) of lower than 80%. In order to validate the cal- 
culations, the simulation results from a previous literature [19] at 
a condition close to our work (Us = 85%) are included in Fig. 2. 
Our calculations show good agreement within a range of high 
operating voltages (0.65-0.75 V). However, at lower voltages, 
the data from the literature shows higher power density. This 
is probably due to the observed temperature increase (within 
100 K) which is particularly pronounced at high current density 
(low voltage) whereas our calculations were based on the isother- 
mal condition. From Fig. 2, regarding the electrical efficiency, 
it is favorable to operate the SOFC at high voltage in order to 
obtain high efficiency. However, when taking into account the 
power density, operation at a high value of voltage is not practi- 
cal due to the achievement of low power density. Therefore, in 
practice suitable operating voltage and fuel utilization should be 
carefully selected. Some workers suggested to operate the cell 
at 70% of maximum power density [20] and fuel utilization of 
80-95% [21-24]. In the present work, the fuel utilization (Uf) of 
80% was considered and SOFC with non-uniform potential oper- 
ation (SOFC-NUP) was investigated with the aim to improve the 
power density without lowering the electrical efficiency. 

Fig. 3 shows the relationship between power density and elec- 
trical efficiency of a simple SOFC-NUP with a cell divided into 
two sections of equal range of fuel utilization (Sp 1 =Sp,2 =0.5). 
Vı and V2 represent the operating voltages of sections 1 and 
2, respectively. The thick solid line indicates the results of the 
SOFC-UP. It is indicated that there are some ranges of opera- 
tion in which the SOFC-NUP offers higher power density than 
the SOFC-UP without lowering the electrical efficiency (area 
above of the thick solid line). The improvement becomes sig- 
nificant when the cell is operated at high electrical efficiency. 


Electrical efficiency [%] 
[z-u9 m] Ajisuap mod 


Electrical efficiency 


Power density 
--Æ&-- Achenbach U,;=85% 


%3 0.6 0.7 0.8 
Operating voltage [V] 


Fig. 2. Performance characteristic curves of typical SOFC-UP (H20/CH4 
ratio =2.2 and T= 1173 K). 
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Fig. 3. Relationship between power density and electrical efficiency of 
SOFC-UP and SOFC-NUP (U;= 80% and T=1173 K; for SOFC-NUP: n=2, 
Sp,1 =0.5). 


In order to indicate suitable operating voltages which offer the 
highest average power density, the value of electrical efficiency 
was specified and the value of Vı was varied. Then V2, which 
gives the desired value of the electrical efficiency, and its cor- 
responding power density can be calculated. For example, as 
shown in Fig. 4, for an electrical efficiency of 43%, the values 
of V2 are 0.623, 0.673 and 0.723 V for V; =0.7, 0.75 and 0.80 V, 
respectively, and the corresponding values of the power density 
are 0.178, 0.156 and 0.159 W cm™?, respectively. By varying 
V,, the maximum power density of 0.162 W cm~? is obtained at 
Vı =0.773 V and V2 =0.700 V. 

According to the above study, the section splits were main- 
tained at S, | =Sp)2=0.5. Those values were then adjusted in 
order to achieve better performance. Fig. 5 shows the effect of 
section split (Sp,1) on the power density improvement at differ- 
ent values of electrical efficiency. It should be noted that the 
reported values are based on the operation using optimum volt- 
ages for each value of the section split. It was found that the 
power density improvement as high as 9.2% can be achieved at 
the electrical efficiency of 45% but the improvement becomes 
less significant when the SOFC is operated at lower electrical 
efficiency. In addition, the optimum Sp; for all cases was found 


48 V,=0.85 ¥~4 0.22 
ee Vi=0.80 y4 0.2 
= 46 seat a 
x 2-77” V,=0.75 V4 0.18 F 
m -- oon = 
is) = ger g 
g 44} =0.80 =- v=0.70 y7 0.16 a 
3) -7 5 
Ea} see 40.14 £ 
D z 
= 42 = 
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E Q 
2 3 
= 40 701 ve 
ore = Electrical efficiency 0.08 
3 Power density 
= 0.06 
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Fig. 4. Effect of operating voltages on performance of SOFC-NUP (n=2, 
Uf =80%, T= 1173 K, Sp. =0.5 and Sp.2 =0.5). 


S. Vivanpatarakij et al. / Journal of Power Sources 167 (2007) 139-144 143 


Table 3 


Comparison of power density between SOFCs with different number of section and section splits 


Electrical efficiency (%) Power density (W cm~?) 


n=1 n=2 (Sp1 =Sp2 =0.5) 


n=2 (Sp,1 =0.55, Sp = 0.45) n=3 (Sp,1 =Sp2 = Sp3 = (1/3)) 


35 0.246 (0.598 V) 
37 0.229 (0.632 V) 
40 0.196 (0.683 V) 
43 0.154 (0.734 V) 
45 0.120 (0.768 V) 


0.248 (0.619, 0.571 V) 
0.232 (0.669, 0.598 V) 
0.201 (0.721, 0.649 V) 
0.162 (0.7726, 0.7001 V) 
0.131 (0.8069, 0.7342 V) 


0.249 (0.631, 0.560 V) 
0.232 (0.666, 0.594 V) 
0.201 (0.717, 0.645 V) 
0.162 (0.768, 0.697 V) 
0.131 (0.803, 0.731 V) 


0.250 (0.644, 0.608, 0.544 V) 
0.233 (0.682, 0.632, 0.583 V) 
0.202 (0.740, 0.666, 0.647 V) 
0.164 (0.782, 0.736, 0.688 V) 
0.134 (0.817, 0.772, 0.720 V) 


10 €=45% 


Power density improvement [%] 
N a 
ee Tae My ae Us 


0 0.2 0.4 0.6 0.8 
Sp, 1 [-] 


Fig. 5. Effect of section split on power density improvement of SOFC-NUP 
(n=2, Up =80% and T= 1173 K). 


to be around 0.55. Table 3 summarizes the values of the power 
density and the corresponding optimum voltages at different 
electrical efficiency of the SOFC-UP and the SOFC-NUPs with 
Sp,1 = Sp,2 and with the optimum Sp,ı and Sp,2. It was found that 
the results of the SOFC-NUP with the optimum Sp; and Sp,2 are 
not significantly different from those of the SOFC-NUP with 
Sp,1 =Sp,2 although the voltages are different. 

To further enhance the performance of the SOFC-NUP, the 
number of separated section (n) was increased. Fig. 6 shows the 
effect of the number of section on the power density improve- 


Power density improvement [%] 


N 


2 3 4 
Number of stage 


wn 


Fig. 6. Effect of number of stage on power density improvement of SOFC-NUP 
(Up = 80%, T= 1173 K and Spx = 1/n). 


ment for three values of electrical efficiency. The section splits 
of all sections are specified to be at the same value of Sp, = 1/n 
in order to simplify the calculations. The optimum voltages of 
the SOFC-NUPs which offer the highest power density for each 
case were determined by following the procedure described ear- 
lier for the case with n=2. It was assumed that the maximum 
power density achieved from the case with adjustable section 
splits does not significantly differ from that achieved from the 
case with equally divided sections. This assumption is valid at 
least for the case of the SOFC-NUP with n=2 as shown ear- 
lier. The calculated results indicate that the obtained maximum 
power density increases with increasing the number of sections 
(n). However, the improvement becomes less significant after 
n>3. The SOFC-NUP with n=3 is likely to be a suitable sys- 
tem for improving the power density without the reduction of 
electrical efficiency. Table 3 also shows the values of the maxi- 
mum power density and the corresponding voltages at different 
values of electrical efficiency. At the electrical efficiency of 45%, 
the power density of the SORFC-NUP with n=3 is 11.7% higher 
than that of the typical SOFC-UP. 

From the above results, it has been demonstrated that the 
use of the non-uniform potential operation with SOFC is techni- 
cally feasible. The cell area can be reduced without lowering the 
electrical efficiency. However, the SOFC-NUP would require 
more sophisticated cell arrangement, power conditioning sys- 
tem and so on. The system control would inevitably become 
more complicated. Further investigations are required before 
implementing this system for commercial use. 


4. Conclusions 


An SOFC-NUP can provide higher power density than a typ- 
ical SOFC-UP without a reduction of electrical efficiency. The 
optimum SOFC-NUP was determined by allowing the operat- 
ing voltage and section split of each section to be appropriately 
adjusted to achieve the highest power density for each level of 
electrical efficiency. The maximum power density can be fur- 
ther improved by increasing the number of separated section 
(n) of the cell; however, it became less pronounced after n > 3. 
Although it is obvious that the non-uniform operation can allow 
the SOFC to be operated at higher performance, further investi- 
gation is necessary to determine whether the cost reduction by 
the reduced stack size would be sufficiently attractive compared 
to the increases of power conditioning cost and complication of 
the SOFC operation. 
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